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Abstract 

The problem of diagonahzation of Hamiltonians of A^-dimensional boson 
systems by means of time-dependent canonical transformations (CT) is con- 
sidered, the case of quadratic Hamiltonians being treated in greater detail. The 
unitary generator of time-dependent CT which can transform any Hamiltonian 
to that of a system of uncoupled stationary oscillators is constructed. The close 
relationship between methods of canonical transformations, time-dependent in- 
tegrals of motion and dynamical symmetry is noted. 

The diagonalization and symplectic properties of the uncertainty matrix 
for 2N canonical observables are studied. It is shown that the normalized un- 
certainty matrix is symplectic for the squeezed multimode Glauber coherent 
states and for the squeezed Fock states with equal photon numbers in each 
mode. The Robertson uncertainty relation for the dispersion matrix of canon- 
ical observables is shown to be minimized in squeezed coherent states only. 



1 Introduction 

The method of canonical transformations (CT) proved to be a fruitful approach in 
treating quantum systems. It is most efficient for systems that are described by 
Hamiltonians, that are quadratic in coordinates and moments, or equivalently in 
boson creation and annihilation operators (quadratic Hamiltonians). The main ad- 
vantage of the method of CT consists in reducing the Hamiltonian H of the treated 
system iS to a Hamiltonian H' of some simple system S' with known solutions. The 
well known example (and probably the first one) of such an application is the di- 
agonalization of the modeled quadratic Hamiltonians in superfiuidity and supercon- 
ductivity theory by means of linear time-independent transformations of boson or 
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fermion operators (the Bogolyubov transforms) [0. In time-dependent CT for 
quadratic systems were used (probably for the first time) in construction of integrals 
of motion that are linear in coordinates and moments. 

Quadratic Hamiltonians model many quantum (and classical) systems: from free 
particle and free electromagnetic field to the waves in nonlinear media, molecular dy- 
namics and gravitational waveguide ^ ^. A considerable attention to quadratic 
classical and/or quantum systems is paid in the literature for a long period of time 
(see, for example, |0, ||, ||, ||, ^ and references therein). 

Diagonalization problem of quadratic Hamiltonians is considered in a number 
of papers 0, |Tl|, [13|, |14|, |^. In general, quadratic Hamiltonians can not be 



diagonalized by means of time-independent CT, even in the one- dimensional case [11 



T3| . In the one-dimensional case the term proportional to the product of coordinate 
and moment can be eliminated by time- dependent CT only. For this purpose a time- 
dependent point transformation (i.e. scale or squeeze transformation) is sufficient 



[TT|. Time- dependent CT are very powerful. Seleznyova has shown that the 
Hamiltonian of a nonstationary quantum oscillator can always be brought to the 
diagonal form of that of the stationary harmonic oscillator by means of linear time- 
dependent CT. 

The aim of the present paper is to establish the canonical equivalence of A^- 
dimensional quantum systems and to perform it explicitly in the case of systems 
with quadratic Hamiltonians. Two systems are called canonically equivalent if their 
Hamiltonians can be related by means of a CT. Due to the known fon Neumann 
theorem CT in quantum mechanics are generated by unitary operators. Therefore 
canonical equivalence is in fact unitary one. A second aim of the present paper is to 
consider the symplectic properties of the uncertainty matrix for canonical observables 
and its diagonalization using linear CT [|l^, 1^ . 



The paper is organized as follows. In section II we show that any two A^- 
dimensional quantum Hamiltonians (time-dependent, in general) H{t) and H'{t) can 
be canonically related via time-dependent unitary operator U{t). The group of CT 
which leave H invariant (i.e. H = H') is shown to coincide with the dynamical 
symmetry group of the system. In the case of two quadratic Hamiltonians the op- 
erator U{t) is an exponent of a quadratic form of coordinates and moments (that 
is, an element the methaplectic group Mp{N, R)). In particular, such operators can 
diagonalize any quadratic Hamiltonian. We note that there are two types of diag- 
onalizations depending of the type of the canonical variables in which the target 
Hamiltonian is diagonal. 

In section III we perform the diagonalization of A^-dimensional quadratic Hamil- 
tonian, expressing the parameters of the corresponding linear CT in terms of solutions 
of linear first order differential equations. For A^ = 1 these equations are reduced to 
the equation z + fl'^{t)z = of classical oscillator with varying frequency. The relation 
of CT to the linear integrals of motion is briefly discussed. 

In section IV the main properties of the uncertainty matrix (T for A^ observables 
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are considered. It is shown that for canonical observables the uncertainty matrix is 
positive definite and thus (due to the known theorem by Wilhamson [|1^, [1^) can be 
diagonahzed by means of hnear CT. For squeezed canonical coherent states (CS) 
and for squeezed Fock states with equal boson/photon numbers in every mode the 
matrix G (when normalized to unity) is found to be symplectic itself. The symplectic 
character of the normalized uncertainty matrix in squeezed CS can also be inferred 
from the results of paper |18 



2 Unitary equivalence of quantum systems 

The main aim in the method of CT is to reduce the Hamiltonian H of the treated 
system 5 to a Hamiltonian H' of some simple system S' with known solutions. CT in 
quantum theory are generated by unitary operators f/, which is called the generator 
of CT. If CT is time-independent then H and H' are unitary equivalent and their 
spectrums are the same. However not any pair H and H' can be related by means of 
time-independent CT. In particular, not any quadratic Hamiltonian can be reduced to 



that of a harmonic oscillator by means of time- independent CT |jTT], |T2|, [T^, |TJ] , even 
in the one- dimensional case |Tl|. The time- dependent CT are much more powerful as 
we shall see below. 

Let |\l/(i:)) be a solution of the Schrodinger equation [ihd / dt — H]\^ {t)) = 0. Then 
for any unitary operator U{t) the transformed state |\l/'(t)), |^l/'(^)) = ^(^)l^(^)); 
obeys the equation [ihd/dt — H']\'$'{t)) = with the new Hamiltonian H', 

H' = U{t)HU\t) - inU{t)dU\t)/dt. (1) 

Conversely, if two Hamiltonians H and H' are related by means of an (unitary) 
operator U{t) in accordance with eq. (P then any solution |\I'(t)) of the system S is 
mapped into a solution \^'{t)) of the system S' . However, not any two given solutions 
|\l/(t)) and \^'{t)) of the two systems could be mapped into each other by means of 
U{t) since U{t) in general cannot act transitively in the Hilbert space. A more 
compact form of relation (D is D\t) = U{t)D{t)U^ (t), where D{t) = ihd/dt - H(t). 
U{t) is intertwining operator for D{t) and D'{t). When D{t) and D'{t) act in the 
same Hilbert space one says that D{t) and D'{t) are unitary equivalent. D{t) is often 
called Schrodinger operator. 

From the requirement for the mean values of the " old" operator A and the " new" 
one A', 

(*(t) |A|^(t)) = (^'(t) |A'|^'(t)), 

it follows that the operators A and A' are related as A' = U{t)AW{t). Therefore the 
new canonical operators of the coordinates and moments and p',^, k = 1, . . . , N are 
related to the old ones as 

q', = U{t)q,U\t), p'^ = U{t)puU\t). (2) 
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Two quantum systems should be called canonically or unitary equivalent if their 
Schrodinger operators are unitary equivalent. The corresponding Hamiltonian op- 
erators H and H', related in accordance with eq. (|lD, should be called canonically 
equivalent with respect to U(t). Let us note the main three advantages of establishing 
unitary equivalence of two systems (see also [||, where in fact canonical equivalence 
of one dimensional oscillators with constant and time- dependent frequencies was con- 
sidered) : 

(a) If we know solutions for one of the two canonically related systems we can 
obtain solutions for the other one as f/(t)|\l/). 

(b) If a time-dependent state |^''(t)) of the system S' is an eigenstate of an operator 
A' then its t/(t)-partner = f/^(t)|^'(t)) in the system S is an eigenstate of the 
operator A = Wlt)A'U{t). 

(c) If the operator A' is an integral of motion for S', i.e. if A' commutes with the 
Schrodinger operator, dA' /dt — {i/h)[A\ H'] = 0, then the operator A = W {t)A'U{t) 
is an integral of motion for the old system S, 



dA/dt-{i/h)[A,H] =0. 



(3) 



This property is very important since if we know one solution for a given system S we 
can construct new solutions acting by the invariant operators on the known solution. 

Proposition 1. Any two N -dimensional quantum Hamiltonians H and H' are 
canonically equivalent. The unitary operator U{t), that relates H and H' takes the 
form 



U{t) = Texp 



to 



H'{t)dt 



UqT exp 



h 



to 



H{t)dt 



= S\t)UoS\t), (4) 



where Uq is constant unitary operator and T and T stand for chronological and an- 
tichronological product. The solution ^ is unique for any initial condition f/(0) = 

Proof Let us perform two successive time-dependent CT by means of t/i = 
UoS^{t) and U2 = S'{t), 



S\t) = texp 



h 



H{t)dt 



(5) 



Then from eq. (|l|) (taking into account dUi{t)/dt = {—i/h)HUi) we easily get 
Hi = for any Uq- The second transformation by means of U2 = S'{t), 



S'it) = Texp 



4 t H'it)dt 
n J 



(6) 



then yields the required result {dUl/dt = dS'^dt = {i/h)S'H'^: 
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H2 = U2H1UI - ihU2dUl/dt = -ihS'dS'ydt = H'. (7) 

Now we see that the direct CT: H ^ H' is performed by the unitary operator (^). 

For a given H and H' the intertwining operator U (t) is not unique. However the 
time-dependence of U{t) is uniquely determined by any initial condition U{0) = Uq. 
Indeed, suppose there is another unitary operator f/(t), which also relates H and H' 
canonically and t/(0) = Uq. Now we note that (it is easily derived from (|T])) if U 
transforms H into H' then U'^ transforms H' back into H and therefor the product 
V — U'^U keeps H invariant: 

H = VHV^ + ^[dV/dt]V^ and V(0) = 1. 

On the other hand, by using eq. (|I]) for U and f/, one obtains the equahty dV/dt — 
(i/^)[V, H] = 0, which means that V is an integral of motion for the system S. Any 
invariant operator for H has the form (note that S{t) is the evolution operator for 
S) V{t) = S{t)V{0)S^{t), and since ^^(O) = 1 we have V{t) = WU = S{t)S^t) = 1. 
In a similar way one can get U{t)U^{t) = 1. And if U'^U = 1 = UU^ , then U = U 
(because of the uniqueness of the inverse U~^). End of the proof. 
Let us note that not necessarily H{0) = H'{0): we have 

H'{0) = UqHUI - inu{<d)[duydt]\t=o- 

Suppose now that H{t) and H'{t) are elements of a Lie algebra C Then S G 
G 3 S', where G is the Lie group generated by C Thus, the CT generator U{t) G G 
(for Uo = 1 and for Uq & G as well) and one can use the known properties of G to 
represent U{t) in other factorized forms. 

The operator (^ converts canonically any A^-dimensional H into any desired 
A^-dimensional H'. In particular H can be converted into H' for a system of free 
particles or for a system of uncoupled harmonic oscillators (A^-mode free boson field). 
In the latter case if if is a quadratic form in terms of A^ canonical operators qk and 
Pj the operator (^) solves the diagonalization problem for quadratic Hamiltonians. 

A CT will be called diagonalizing if the new Hamiltonian H' in terms of the co- 
ordinates and moments is diagonal quadratic form with constant coefficients, i.e. H' 
is a Hamiltonian for a system of uncoupled harmonic oscillators i^ho- One has to 
distinguish between two different kinds of diagonalization of H\ 

First kind diagonalization: H' is diagonal in terms of the new variables q'jiP'k-, 
Second kind diagonalization: H' is diagonal in terms of the old variables qj,Pk- 

In the first case the two systems S and S' are treated in two different (g- and q'- ) co- 
ordinate representations (wave functions "^{qjt) = (g|^'(t)) and \E''(g',t) = (g'|\E''(t))), 
whereas in the second case one can work in the same g-representation (wave functions 
^{q,t) = {q\^{t)) and ^'{q,t) = {ql^'it))). 
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The second kind diagonalization is achieved by means of operator U{t), eq. @, 
with H' of the form of Hamiltonian of N uncoupled stationary oscillators (in terms 
of old variables), 



1 



N r 



2 , 



1 



TUk 



H' = -Y. —vl + ^ku^hl =HUp,q), (8) 



The target Hamiltonian H' may also be taken as a sum of stationary oscillators ifho 
in terms of the intermediate variables ql^\ p^j^^ as well. In the latter case the second 
CT {£\ p'i^) {qip',), generated by U^it) = exp[-{z/h)HUq^'\p(^y)t] = S^t), 
takes the explicit form of rotations 

Qk = Qk^ cos(^fct) + ^P^^ sin(^fct), 

Let us briefly elucidate the two CT involved into the proposition 1. The first 
one, generated by Ui = UoS^{t), brings H to zero, therefore the new states |\Ef)i are 
time-independent. This is because S^(t) is an evolution operator for the S back- 
ward in time. After the first CT (generated by Ui) the new canonical variables 
g(^) = Ui{t)qkUl{t) and pf^ = Ui{t)pjUl{t) obey the equations {dUi/dt = iUiH, 
dUl/dt = -iHUl) 

^ = i|c/.i/f/;.,<"i. ^ = i|c/,Hf/f.p«], (10) 

i.e., q^^\ p^^^ are Heisenberg operators for the old system S. 

The generator of the second CT U2{t) = S'(t) is recognized as the evolution 
operator forward in time for the target system S'. In the construction (^ U2 is 
applied to the intermediate Hamiltonian Hi. 

It is worth noting at the point the case of CT in the system S, generated by 
its own evolution operator S{t). This CT converts H{t) into Hamiltonian H"{t) = 
S{t)H{t)S^{t)+H{t). If H is time-independent then S{t)HS\t) = H and H" = 2H. 
From (^'(t)|A|^(t)) = (^'|5''f(i:)AS'(t)|^') we derive that the new canonical variables 
in this case, 

q'l = S{t)qkS\t) ^ ql pi = Sit)p,S\t) ^ pi (11) 

when expressed in terms of the old ones, g^, pj, are integrals of motion of S, satis- 
fying the eq. (§). Such integrals of motions for quadratic systems H{t) have been 
constructed in and intensively used later p, ||, |[ . 

Consider the symmetry of H under CT. We want to specify the set of CT for 
which H', defined in (^, coincides with H, i.e. we look for CT that keep H invariant 
(and thus keep the Schrodinger equation invariant), 

H' = U{t)HU\t) - ihU{t)^^^^ = H. (12) 
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For time- independent U eq. (^) reduces to H = UHW . From (0) For H' = H the 
CT generator is (see (H)) U{t) = S(t)UoS'^ (t), where Uq is arbitrary unitary operator. 
Then dW{t)/dt = {i/h)[W (t), H] and we see that the equahty in ( |T2|) is identically 
satisfied. Thus, the CT generators U{t) for which H' = H have the form S'(t)f/o5'"^(t), 
i.e. U(t) are integrals of motion for the system: [U{t),D(t)] = 0, where D{t) is the 
Schrodinger operator, D{t) = ihd/dt — H. In the first paper of refs. [0 the dynamical 
symmetry group of a system S has been defined as a group of unitary operators, that 
commute with D{t) and act irreducibly in the Hilbert space. Now we see that this 
symmetry group leaves H' = H and is highly nonunique, since the unitary operator 
Uq in U{t) is arbitrary - one can take Uq from irreducible representations of any Lie 
group. Then the set of invariants S(t)UoS'^ (t) realize an equivalent representation of 
the same group. For example, by means of the invariants g° and p° one can construct 
an irreducible representation of the Lie algebra of the Heisenberg-Weyl group H\y{N) 
and the quasi unitary group SU{N, 1) as well 0. This means that the groups H\y{N) 
and SU {N, 1) can be considered on equal as dynamical symmetry groups of any A^- 
dimensional system. 

In the next section we consider the above described unitary (canonical) equiv- 
alence approach in greater detail for quadratic quantum systems, for which some 
explicit solutions can be obtained. 

3 Canonical transformations of quadratic systems 
and diagonalization. 

We consider the general A^-dimensional nonstationary quantum system with Hamil- 
tonian H{t), that is a homogeneous quadratic form of coordinates and moments, 

H(t) = Ajk(t)pjPk + Bjk(t)pjqk + Bjk{t)qjPk + Cjk(t)qjqk, (13) 

where the coefficients Ajk(t) = Akjit), Bjk(t), Bjk(t) and Cjkit) = Ckjit) are arbitrary 
functions of time. From W = H it follows that Ajk{t) and Cjk{t) are real, and 
Bjk(t) = Bljit). It is not a significant restriction to take Bjk real and put Bjkit) = 
Bkj (t) (the imaginary parts of Bjk can be eliminated by adding a non-operator term to 
H). In (|T3p the summation over the repeated indices is adopted. We can introduce 
A^-component vectors q = (gi, q2, ■ ■ ■ , g^v), p = {pi,P2, ■ ■ ■ ,Pn), and N x N real 
matrices A(t), B(t), C{t) {A(t) and C(t) are symmetric) and rewrite the Hamiltonian 
(|l^) in a more compact form 

H{t) = pA{t)p + pB{t)q + qB{tfp + qC{t)q, 

where is the transposed of B. To shorthand the notations it is convenient to 
introduce the 2iV-vector Q = {p,q) and 2N x 2N matrix Ti. (the grand matrix) and 
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rewrite the Hamiltonian (0) as {fijU = 1,2,..., 2N) 



H{t) = n,,,Q^Q. = Qn{t)Q, n=[ n ] ■ (14) 




We note that nonhomogeneous quadratic Hamiltonians (i.e., Hamiltonians of the 
form (0), (|I4|) with hnear terms added) can be easily reduced to the forms (|13|), (0) 
by means of simple time-dependent displacement transformations. 

Let H' be an other quadratic Hamiltonian 



H'it) = Qn\t)Q, n' = c' )■ 

Then the unitary operator U{t), eq. (^, which relates canonically Hamiltonians (p^) 



and ([T5|) , is an exponent of a quadratic in q and p form (we take Uq G Mp{N, R)), 

U{t) = S'{t)UoS\t) = exp 



(16) 



where ^{{t) is a new grand matrix of the form (|1^) and (|T^. 7^(t) can be expressed 
in terms of the Hamiltonian matrices ?^(t) and Ti^t) using the Baker- Campbell- 
Hausdorff formula. In this case the operator ( |16|) generates linear transformation of 
coordinates and moments (we write it in x and 2N x 2N matrix forms), 

«' - ( r ) - ( t t ) ( ? ) • 

where Xpp, Xpq, Xqp and Xqq are N x N submatrices of A(t). 

From eqs. ([l|), (p!^ , (|15]) and (|T7|) we obtain the following relation between the 



symmetric matrices H, Ti' and 7i (0) and the symplectic matrix A, 

jH{t) = -H'{t) + A^n{t)A. (18) 

We see that for a given 7i{t) and 7i(t) this is a simple linear equation for T-C'{t). 
However for a given Hamiltonian matrices Ti. and Ti' this is highly nonlinear equation 
for 7i(t) since the matrix A{t) is to be expressed in terms of 7i(t) again: AQ = 
U{t)QU\t). Nevertheless for any given (differentiable with respect to t) matrices 
7i(t) and 7^'(t)) and for a given initial condition Ti^ the above system of equations has 
unique solution for ?i(t), since the expression of A in terms of Ti is also differentiable 
and Peano theorem could be apphed . 

In this scheme A(t) is naturally represented as a product of two other 2N x 
2N matrices A*^^) and A^^-* of the form ([T7| ) corresponding to the two successive CT 
generated by Uiit) and U2{t): 

A = A(2)A(i); g(i) = A(i)g, Q' = A^^^Q^^\ (19) 
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The matrices A*^^-* and A^^-* are seen to be solutions of the first order hnear equations, 
d^d) ^ A(i)ir(i)(t), ^A(2) = F(2)(t)A(2), (20) 

where 

F«(t) = -2J7^(^), F^'\t) = 2jn'it), -^=[-1 J ) • (21) 

If H' is diagonal as for the oscillator system (|^) then the second eq. (|20D is easily 
solved: A*^^^(t) = exp(2J7iho^)Ao^''. To perform the diagonalization of a quadratic H 
one has also to solve the first matrix equation in (pO]) and obtain A^^^ (t), which in prin- 
ciple is always possible. In the case of stationary initial H the T exponent becomes or- 
dinary one, so the explicit solution is given by the matrix exponent Aq^"* exp(— 2J7it). 
So for stationary H the total A matrix takes the form 

A{t) = exp{2Jnhot)Ao^Ao^ exp(-2J7^^), (22) 

where Aq-* are arbitrary symplectic matrices. One can put Ag^'^^ = 1, which corre- 
sponds to ?7o = 1 in eq. (H). Having obtained explicitly A(t) one can next try to 
solve eq. (|TB|) and obtain the generating operator U{t) in the form of the quadratic 
exponent (piBI) . 

Note, the resulting H' is diagonal in the variables, which we choose for ifho- Let 
those variables be , ql^\ Then the final variables p^, obey eqs. (^. Inverting 
the transformations (^ we obtain H' diagonal in terms of the final variables as well: 
H' = Hho{p', q'). In this way we perform explicitly the first kind diagonalization. If 
H' = ifho in terms of old variables pk, qk (second kind diagonalization), then H' is 
evidently not diagonal in terms of p^, g^. 

For some time-dependent H{t) explicit solutions of eqs. (pOD can also be found. 
Thus, in the case of = 1, following the scheme of refs. 0, |^, one can express matrix 
elements of A^^\t) in terms of a complex function z{t)^ that obeys the equation of 
classical oscillator z^-VL^{t)z = 0, where f2^(t) is simply determined by the parameters 
A, B, C of the Hamiltonian (|l^) (for = 1 these are not matrices, therefore we put 
A = a, B = b, C = c), 

fi2(t) = 4ac + 2ba/a + a/2a - 3a'^/Aa'^ - 46^ - 2b. 

For harmonic oscillator with varying frequency uj{t) we have Q'^{t) = uj^{t). It is 
seen that an Q{t) corresponds to a class of quadratic H{t). For example constant 
Q corresponds to the stationary oscillator and to the oscillators with varying mass 
(damped oscillators) m(t) = niQ exp(— 26t) and m(t) = ttiq cos^ bt, considered later by 
many authors (see refs. in ^ ^). Analytical solutions to the equation of z{t) are 
known for a variety of "frequencies" Q{t). In the case of an oscillator with varying 
frequency the diagonalizing CT generator U{t) has been expressed in terms of z{t) 
in §. 
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Let us briefly discuss the algebraic properties of the matrix A{t) and its subma- 
trices Xpp, Xpg, Xgp, and Xqq. From the canonical commutation relations it follows that 
A(t) obeys the relation (the symplectic conditions, J defined in eq. (|2l|) ) 

AJA^ = J, (23) 

which for the N x N matrices Xgg, Xpp, Xgp and Xpg, defined in eq. ([T7|)), read 

\ \T _ \ \T _ -I \ \T _ \ \T \ \T _ \ \T /r)A\ 

^PP^qq '^Pq^qp ^i ^q<i^qp ^QP^qqi ^Pl^pp "^PP^pq- l^^/ 

The set of matrices that obey the relation (^) is defined as the symplectic matrix 
group Sp{N, R) (the transformation x' = A'^x preserves the quadratic form xJx). It 
has A^(2A^+ 1) real parameters. The rank of its Lie algebra is (following |jl9| we use 
the notation Sp{N, R) instead of Sp{2N, R)). It is known that in classical mechanics 
the set of linear homogeneous CT forms a symplectic group Sp{N, R) . In the quantum 
case the set of matrices A, that realize homogeneous linear transformations of the 
operators of coordinates and moments close the same group. However the set of 
unitary operators U for which UqW and UpW are linear combinations of p and q 
contains one extra parameter, namely the phase factor. If one considers CT in greater 
detail as transformations of coordinates, moments and vectors in Hilbert space one has 
to count the phase factors as well and then we get the larger group Sp{N, R) x U{1) = 
Mp{N, R). If we consider transformations of coordinates, moments and states we have 
to factorize over [/(I): Mp{N, R)/U{1) = Mp{N, R). The resulting group Mp{N, R) 
is called methaplectic group. It is double covering of Sp{N,R). The Lie algebras of 
Mp{N, R) and Sp{N, R) are isomorphic [|1^, |2^. They are of dimensions N{2N + 1) 



and this is the number of independent matrix elements of matrix Ti. in (|16]). The 
generators U(t) of linear CT (|T^ can be considered as operators of the unitary 
(but not faithful) representation U{A) of the symplectic group Sp{N,R). One can 
use the group representation technique [|l9l to represent U(t) G Sp{N, R) in several 



factorized forms. In the case of one dimensional nonstationary harmonic oscillator 
the diagonalizing CT generator U{t), U{t) G SU{1, 1), and its factorized forms have 
been considered in [Rl. 



If one considers Hamiltonians ([T3|) with linear terms d{t)p + e{t)q added, then in 
the same way one would get that such inhomogeneous quadratic Hamiltonians can 
be diagonalized to the form (H) by means of the same U{t), eq. (^, this time U{t) 
being an element of the semidirect product group Mp{N, R)k> if^(A^), where if^(A^) 
is the A^ dimensional Heisenberg-Weyl group. 



4 Diagonalization of uncertainty matrix and min- 
imization of characteristic inequalities 

The established possibility of converting (by means of time-dependent CT) any A^- 
dimensional Hamiltonian H to that of the system of uncoupled harmonic oscillators 
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suggests to expect that the dispersion matrix (j{Q,p) of canonical observables , 
V = 1, . . . 2iV, in any (generally mixed) quantum state p could be diagonalized by 



means of some state dependent CT. It turns out that this really holds |jT5|, |T6|. 

Let us recall the notion of dispersion matrix (7(X,p) (called also fluctuation ma- 
trix, or uncertainty matrix). This is an x n matrix constructed by means of the 
second moments (the variances and covariances) of observables Xi, . . . , X„ in a state 
p. The matrix elements (J^i, of (7 are defined as covariances AX^j^Xy of the observables 
X^ and Xy, i^, p = I, . . . ,n, 

a,,{X,p) = i(X^X, + X,X^) - (X^)(X,) = AX^X,(p). 

The matrix (j{X,p) is symmetric by construction. It satisfy the characteristic uncer- 
tainty relations ^ G("H(t(X,p)) > C('")(C(X,p)), where C(X,p) is then anti- 
symmetric matrix of the means of commutators of X^ and X^,, C^^ = — z([X^, X;^])/2, 
and C^")(M), r = 1, . . . ,n, are the characteristic coefficients of a n x n matrix M 
2^ . The characteristic coefficient of maximal order r = n is the determinant of M. 



The characteristic uncertainty relation of maximal order r = n, 

det cr(X, p) > det C(X, p), (25) 



has been established by Robertson and is called Robertson uncertainty rela- 
tion. For N = 2 inequality ( P^D recovers the Schrodinger uncertainty relation 
(AX)2(AF)2 - (AXr)2 > |([X, Y])f /4, which for the canonical pair [q,p] = i, 
takes the simpler form of (hereafter we put h = 1) 

{Apf{Aqf - {Apqf > 1/4. (26) 

The proof of ( ^5]) is based on the nonnegativity of the matrix R = (J + iC [ p]] . 
Properties of R (to be called Robertson matrix) are reviewed in [Q. Here we need 
the nonnegativity property of (7(X,p). 

Proposition 2. The uncertainty matrix for any n observables Xi, . . . , X„ is non- 
negative definite, (7(X,p) > 0. 

Proof. The proof relies to the Robertson inequality (^) and on the observation 
that the principal submatrices m(Xjj, . . . , Xj^, p), r < n, of (7 can be regarded as 
uncertainty matrices for r observables Xj^ , • • • , Xj^ in the same state p. Therefore 
the submatrices m(Xjj, . . . , Xj^), p) also satisfy Robertson relation (p5|), i.e. their 
determinants (the principal minors of (j) are nonnegative. And if all principal minors 



of a matrix M are nonnegative, then M > End of the proof. 



The uncertainty matrix (j{Q,p) for canonical observables Q^, p = 1,...,2N: 
Qk = Pk, Qn+k = Qk, k = ^, ■ ■ ■ , N, possess some further properties. 
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Proposition 3. The uncertainty matrix for 2N canonical observables is posi- 
tive definite, (j{Q,p) > 0. 

Proof. From the canonical commutation relations [qk,Pj] = i^kj it follows that 
det C(g,p) = (1/4)^. Then (||) yields 

det(j(Q,p)>(l/4)^. (27) 

As a symmetric matrix (j{Q,p) can be diagonalized by an orthogonal transforma- 
tion ^ Q'^ = l^uQu- The uncertainty matrix for new observables Q'^ is (j' = 

(7{Q',p) = ■yCTj'^. This transformation preserves the determinant (and all the other 
characteristic coefficients) of (J. In view of det (j{Q,p) > all diagonal elements of 
a' are positive. Therefore (7{Q,p) > 0. End of the proof. 

The desired diagonalization of (j{Q,p) using linear canonical transformations 



Tq , |T6| now follows from the Proposition 3 and the known theorem [0, [T2|, |T3| 
that any positive definite symmetric matrix M can be diagonalized by means of con- 
gruent transformation with a symplectic matrix A, M — M' = AMA^. In [15| 



the diagonalization of (j{Q,p) is performed explicitly by means of three consecu- 
tive linear canonical transformations. The diagonal elements of the diagonalized 
a' are variances (AQ^)^. Additional scaling transformations — > gf = a^g- with 
ai = (Ap-/Agj')^/^ equalize the variances of q" and p" = Pi/ai. Note that: (a) the 
diagonalizing symplectic matrix A is not unique [O, ES; (b) A is state-dependent. 



Therefore it may depend on time when the state is time-dependent. 

Denoting the generator of the total diagonalizing canonical transformation by 
f/(A) [Q' = f/(A)Qf/t(A), f/(A) e Mp{N,R)] we obtain the equality a{Q,p') = 
<t{Q',p), where p' = UpW. Thus every state p is unitary and methaplectically 
equivalent to a state p', in which the uncertainty matrix (j{Q,p') is diagonal with 
equal variances of coordinates qi and moments pf. Aqi = Api. If the initial state 
|\E') is pure time-dependent state of system S with Hamiltonian H, then the CT 
is time-dependent and the new state 1^^') obey the Schrodinger equation with new 
Hamiltonian (|1]). 

Examples of pure states with diagonal uncertainty matrix with equal variances of 
coordinates and moments are Glauber multimode coherent states \a) and multimode 
Fock states |n). Therefore in the Klauder-Perelomov Mp{2,R) CS \g,d) = U{g)\d) 

— * 

and \g,n) = U{g)\n) {g being the group element) the dispersion matrices (j{Q',g,d) 
and (7{Q',g,n) are diagonal and with equal variances of g,' = {g)qiU{g) and p[ = 
U\g)piU{g). In \d) all variances are equal and minimal, Aqi = Api = l/V^, whereas 
in Fock states the variances are equal in pairs, (Agj)^ = {ApiY = 1/2+nj. Multimode 
CS I a) minimize Robertson inequality (P7|), whereas in \n) one has det (T{Q,n) = 
n.(l/2 + ni). 

It is clear from the above consideration that the uncertainty matrix in any group- 
related CS T(5f)|\I'o) with reference vector |\Efo) equal to \d) or \n) is diagonalized 
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by the CT Q'^ = T^{g)Q^T{g), that is hnear for the group Mp{N,R) only. In 
physical literature the Mp{N, R) group-related CS U{g)\d) and U{g)\n) are known as 
multimode squeezed CS and squeezed Fock states respectively. The operator U (g) e 
Mp{N,R) can be called multimode squeeze operator its canonical form being 
exp[(a ''^2:0 ^ — az*a)/2], where a^za'^ = ajzija^j, i,j = 1,...,N p6|. It is more 



adequate to call it squeeze and correlation operator since, e.g., for pure imaginary zu 
it generates covariances of pi and and doesn't squeeze, while for real zu it generates 
squeezing and doesn't correlate. The wave function {x\U{g)\d) of Mp{N,R) CS is 
Gaussian (an exponent of A^-dimensional quadratic form), thereby that states are 
also called Gaussian pure states 



It is interesting to note that the multimode squeezed states U{g)\d) are the unique 
states to minimize the Robertson inequality 



Proposition 4. The equality in the multimode Robertson uncertainty relation, 
eq.(27), holds in the multimode squeezed states U{g)\d) (g G Mp{N, R)) only. 

Proof. Let A(p) be a symplectic matrix that diagonalizes the dispersion matrix 
Cf{Q,p), and U = U{A) - the generator of the diagonalizing CT Q' = A{p)Q = 
uIa)QW{A). U{A) belongs to Mp{N,R). We have 

a{Q', p) = A{p)a{Q, p)A^(p) = a{Q, p'), p' = U{A)pU\A), (28) 

a{Q, p) = diag{si, S2, . . . , S2n}, (29) 

where the diagonal elements s^, are the variances of qk and pk in the new state p': 
Sk = (Apfc(p'))^, SN+k{p') = {^Qk{p')y- The determinant of (j{Q,p') is a product of 
all diagonal elements s^, v = 1, . . . , 2N , 

det Cr((5,p) = detCr((5,p') = [siSAr+i][s2SAr+2] . . . [sArS2Ar]. (30) 

From Heisenberg uncertainty relation we have for every factor in eq. (^) the in- 
equality 

suSN+k = {^Pkfi^qkf > 1/4. (31) 

From eqs. (|3T|) and (|30|) we derive that the equahty in Robertson relation (pT]) holds 
iff the equality in eq. (|3T| ) holds for all modes (for every k = 1, . . . N). The minimal 
value of 1/4 of the product of variances of q and p cannot be reached in mixed state 
- it is reached in the Stoler states |Q \a, r) = exp[r(a^^ — a^)] \a) only (see proof in 
the Appendix). Thus the equality in (p7|) holds in pure states U{A) Yik S{rf:)\d) only. 
The unitary operator S{r) = exp[r(a''"^ — a^)/2] (the squeeze operator) belongs to 
U{1,1) ~ Mp{l,R). Therefore the unitary operator [/(A) Ofc 5'(rfc) = U{g) belongs 
to Mp{N, R), and the unique minimizing states are Mp{N, i?)-group related CS with 
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reference vector \a). End of the proof. 

Since Glauber CS \a) are eigenstates of every annihilation operator (with 
eigenvalues Ofc, = 1, . . . , A^), the minimizing states U{g)\a) are eigenstates of the 
canonically transformed annihilation operators = U{g)akU^{g), which are linear 
combinations of ai, . . . ,ai\f: a'^ = Ukjaj + VkjOj- Therefore the minimizing states 
U{g)\d) (the multimode squeezed states) can be denoted equivalently as \d,u,v). 
For V = (and u = 1) they coincide with \d). 

For quadratic Hamiltonians the time evolution operator Uqua.d(t) G Mp{N,R). 
Therefore the time evolution of \d,u,v) for quadratic Hamiltonians is stable, i.e., 
U{t)\d,u,v) = \d,u(t),v(t)). The evolved states \d,u(t),v(t)) are eigenstates of the 
new annihilation operators Ak{t) = U {t)a'^U\t) , which are again linear in aj and a] 
and are integrals of motion of quadratic system. Overcomplete system of eigenstates 
\d,t) of integrals of motion Ak{t) has been constructed in ref. [5] and used later in 
many papers [10]. 

A further property of the uncertainty matrix (the fourth one) we want to note 
here is referred to its symplectic character, the normalized uncertainty matrix O" = 
(7 / (deter) ^'^^^ is symplectic for a certain class of states. In order to find out that 
states we note the invariance of the symplectic property of a matrix M under the 
congruent transformation AMA'^ with a symplectic A: if M is symplectic, that is 
MJM'^ = J, then M' = AMA^ is also symplectic. This symplectic invariance 
can be easily proved using the known property that if AJA^ = J then one also 
has A^JA = J: M'JM'^ = AMA^JAM^A^ = AMJM^A^ = AJA^ = J. This 
invariance enables us to study the symplectic properties of (7 in its simpler diagonal 
form. For diagonal uncertainty matrix a = diag{si, . . . , sat} the symplectic condition 
(7 J (J = J reduces to 

siSn+i = S2SN+2 = ■■■ = SNS2N = (det cr)^/^ . (32) 



One solution to (|32D can be immediately pointed out, recalling the meaning of 
as the variance of Q^: the uncertainty matrix in the multimode Glauber CS \d) is 
diagonal with Sk = {^PkY = 1/2, s^v+fc = i^Pk)^ = 1/2, k = 1, . . . ,N, which clearly 
satisfy (0). Therefore the normalized uncertainty matrix in pure states U{g)\d) 
that are unitary equivalent to Glauber CS with U{g) G Mp{N, R) is symplectic. 
These states, as we have already noted, are called Gaussian pure states or multimode 
squeezed CS. In fact the symplectic character of the normalized uncertainty matrix 
for Gaussian pure states was established in [|18|: in that states our 5" is equal to 2(7 
and this quantity coincides with the matrix G(U_~^, ~V_) of [jl8|, which was shown to 
be symplectic [ff 



A second solution to ( |5^ is provided by the uncertainty matrix in (multimode) 
Fock states \n) with equal numbers Uk = n (equal numbers of photons in every 
mode). In \n) we have (Ap^)^ = 1/2 + Uk = (Ag^)^- Therefore in states U{g)\n) 
with rii = . . . = TiN and U{g) G Mp{N,R) the normalized uncertainty matrix is 
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symplectic. The above two families of states do not exhaust the set states with 
symplectic (normahzed) uncertainty matrix. 

Let us write down the symplectic conditions and the Robertson relation for 
(j{Q,p) in terms of the four N x N blocks app{p), (Jgq{p), crpg{p) and (Jgp{p), 

^(Q,p)=f"^f| "^fM, (33) 

Inserting this into (7 J (7 = J, and taking into account that app and cr^g are sym- 
metric, and (Tpg = ajp we obtain 

o-ppO-gq - {(Jpqf = (det cr)^/^ , (34) 

^pp^qp ^pq^pp 0) ^qp^qq ^qq^pq 0- (35) 

Squeezed CS U{g)\d) minimize (P7|), i.e. det cr = (1/4)^. Therefore in U{g)\a) the 
symplectic condition (Q) reads (TppUgg — (upg)^ = 1/4. The latter formula was obtained 
in 1^ for the squeezed CS of the form exp[(a^2;a^ — a2;*a)/2] \a) by direct calculations 
(but with no reference to Robertson inequality, neither to the symplecticity of the 
uncertainty matrix). In squeezed Fock states U{g)\n) we have det (T{g, n) = nfe(l/2 + 
^k) > (1/4)^- For these states the symplectic condition (|34D is valid iff = n, and 
reads appagg - {apgf = (1/2 + nf. 



In terms of the N x N matrices Robertson inequality (27) takes the form (using 
known formulas for the block matrices [p^ ) 



det [cTppCrqg - (TppfrgpfTppCrpg] > (1/4)^. (36) 

For a symplectic we have (JppCTgp = Cpgapp, and the Robertson relation simplifies 
to det[crpp(jgg — (cTpq)^] > (1/4)^. This form is quite similar to that of Schrodinger 
inequality (pGf ) for p and q: for = 1 we have apg = Apq, (Xpp = App = {ApY, and 
agg = Aqq = (Aq)^ . 

It is curious to note that the Robertson matrix R for normalized d'{Q,p) and 
C{Q,p) = C{Q,p)/ (det C{Q,p)^ , R = a + iC, is also symplectic for squeezed 
CS and squeezed Fock states with Uk = n: RJR^ = J, that is i? G Sp{N, C). 



5 Appendix 

Proposition Al. Heisenberg inequality (Ag)^(Ap)^ > 1/4 is minimized in the Staler 
states \a,r) = exp[r(a"''^ — a^)/2] \a) only. 

Let p be a general mixed state. Any mixed state can be represented in the form 
P = J2k Pk\4'k){4'k\, where pk > 0, and {I'i/'fc)} is some complete orthonormal set of 
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pure states. The mean value of an operator X in p is given by (X) = Tr(Xp). 
Consider the mean value of the operator 6"''(A)6(A), where 

b{X) = Xq + ip-{X{q) + i{p), XeR. (37) 

For positive A the operator 6(A) is, up to a factor l/v^2A, a boson annihilation 
operator, [6, 6^"] = 2A, and 6^6 is, up to a factor 1/2A, the number operator, which is 
nonnegative definite. For negative A the operators 6 and 6"'" are interchanged. The 
mean of 6"''(A)6(A) in p reads 

(6t(A)6(A)) = ^p,(^fc|6t(A)6(A)|^,), (38) 

k 

where all means {ipklb^ {X)b{X)\i/jk) are nonnegative. On the other hand, by the use of 
(^), this nonnegative mean (6'''(A)6(A)) can be written as 

(6t(A)6(A)) = A'(Ag)2 - A + {Apf > 0. (39) 

The A-roots of the equation A^(Ag)^ — A + (Ap)^ = must be real, wherefrom one 
deduces Heisenberg inequality. The equality in Heisenberg relation corresponds to 
the equality in (P^), i.e. to the vanishing {6^(A)6(A)). From ( pS] ) it is seen that 
(6"''(A)6(A)) = if and only if {ilJk\b^{X)b{X)\ilJk) = for every k = 1,... (in view 
of pk > 0). From the uniqueness of the vacuum state it follows that all pk but one 
(say pi) must be zero. Therefore (6^'(A)6(A)) = in pure state \ip) only and iff it 
is an eigenstate of Ag + ip. The final step is to identify the minimizing pure state 
with \a,r). The minimizing state must be eigenstate of Xq + ip for some real A. 
\a,r) are eigenstates of Xq + ip) with A = [cosh(2r) — sinh(2r)]/[cosh(2r) + sinh(2r)] 
and eigenvalue av^/[cosh(2r) + sinh(2r)]. Thus (Ag)^(Ap)^ = 1/4 holds in states 
|a, r) only. In slightly different notations the proof of the statement that a state with 
absolute minimum of the product (Ag)^(Ap)^ is a pure state is given in [P^]. 
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